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ABSTRACT 

We investigate the most plausible stellar Initial Mass Function (IMF) and the main origin 
of the tilt of the Fundamental Plane (FP) for old, massive early-type galaxies. We consider 
a sample of 13 bright galaxies of the Coma cluster and combine our results with those ob- 
tained from a sample of 57 lens galaxies in the same luminous mass range. We estimate the 
luminous mass and stellar mass-to-light ratio values of the sample galaxies by fitting their 
SDSS multi-band photometry with composite stellar population models computed with dif- 
ferent dust-free, solar-metallicity templates and IMFs. We compare these measurements and 
those derived from two-component orbit-based dynamical modelling. The photometric and 
dynamical luminous mass estimates of the galaxies in our sample are consistent, within the 
errors, if a Salpeter IMF is adopted. On the contrary, with a Kroupa or Chabrier IMF the two 
luminous mass diagnostics differ at a more than 4cr level. For the massive Coma galaxies, their 
stellar mass-to-light ratio scales with luminous mass as the corresponding effective quantities 
are observed to scale on the FP. This indicates that the tilt of the FP is primarily caused by stel- 
lar population properties. We conclude that old, massive lens and non-lens early-type galaxies 
obey the same luminous and dynamical scaling relations, favour a Salpeter IMF, and suggest 
a stellar population origin for the tilt of the FP. The validity of these results for samples of 
early-type galaxies with different age and mass properties still remains to be tested. 
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1 INTRODUCTION 

It is now commonly accepted that the most massive (M, > 
10" Mq) galaxies in the local Universe are mainly ellipticals, com- 
posed of uniform old stellar populations that are completely formed 
at z ~ 1. For these objects, a possible fo rmation pictu re consists in 
a very rapid "monolithic" collapse (e.g.. ^ Larsorjl 1 974 ^ followed by 
some additional moderate star formation and merging activit y, as 
obser ved in ellipticals from z ~ to z ~ 1 (e.g.. lConselice et al.l 
l2007h . Although mass seems to drive for the most p art the galaxy 
star formation history (SFH, e.g.. lCowie et al.iri996l) . environmen- 
tal effects have also observable consequences (e.g., iGobat et"ail 
[2OO8). Understanding the initial formation and the subsequent evo- 
lution of massive galaxies is thus an important test for the ACDM 
paradigm. 

Strong gravitationa l lensing is an invaluable t ool to study the 
internal structure (e.g., Rusin & Kochanek 2005) and the lumi- 
nous and dark matter composition (e.g., Grillo et al. 20 08b. 2009bi) 
of massive early-type galaxies. Several studies (e.g., lOrillo et al.l 
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l2009al . hereafter G09a: IXreu et al.ll2006l) have proved that early- 
type lens galaxies are an unbiased sub-sample of massive early-type 
galaxies. Hence, massive non-lens as well as lens early-type galax- 
ies can be used to conclude on their physical properties, like the 
stellar IMF, and the origin of the observed scaling relations. In par- 
ticular, the FP (Djorgovski & Davies 1987; Dressier et al. 1987), a 
scaling law between effective radius, central stellar velocity disper- 
sion, and average surface brightness, and its interpretation (known 
as the tilt of the FP) in terms of a systematic increase of the galaxy 
effective mass-to-li ght ratio with effective mass have been exten- 
sively studied (e.g.. lCiotti. Lanzoni & Renzin i 1996). However, for 
old, massive early-type galaxies, it is still no t clear which IMF is the 
most suitable one ("e.g.. Ivan Dokkunj|2008l and references therein) 
and if the tilt of the FP is mainly dete rmined by variations in the 
galaxy stellar popula tion features (e.g., Allanson et lj]|2009l) . dark 
matt er content (e.g., 



Tortora et al. 



(e.g.. lBertin. Ciotti & Del Principell 20021) 



200^) 7 or structural properties 



This letter is organised as follows. In Sect. 2, we describe how 
and why we select a sample of massive early-type galaxies of the 
Coma cluster. In Sect. 3, we compare photometric and dynamical 
modelling to measure the amount of luminous mass in the galax- 
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Table 1. Dynamical and photometric luminous mass and stellar mass-to- 
light ratio estimates for a sample of 13 early-type galaxies of the Coma 
cluster. 
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ies and analyse the consequences of adopting different IMFs. In 
Sect. 4, we explore for the galaxies in the sample the relation be- 
tween stellar mass-to-light ratio and luminous mass and discuss 
the implications of this relation in the tilt of the FP. In Sect. 5, 
we combine the results on massive non-lens (Coma) early-type 
galaxies with those obtained in previous works on massive lens 
(SLACS) early-type galaxies. In Sect. 6, we summarise our anal- 
ysis. In this study, a concordance cosmology (Ho = 70 km s"' 
Mpc"', flm = 0.3, Aa = 0.7) is assumed and all logarithms have 
base 10 and adimensional arguments, obtained by dividing the di- 
mensional quantities by the corresponding measurement units. 



2 THE SAMPLE 

We concentrate on a sample of 13 bright early-type galaxies of the 
Coma cluster. The sample is composed of flattened, rotating and 
non-rotating, elliptical and lenticular galaxies for which detailed 
axisymmetric, orbit-based, dynamical modelling of kinemat ic and 
photo metric observations was performed by Thomas et al. l l2007l : 
l2009l hereafter T07 and T09). Among the galaxies studied in T09, 
we select those that have multi-wavelength (iigriz bands) photomet- 
ric measurements from the SDS^ 

The choice of this sample is dictated by several different rea- 
sons. First, we aim at comparing published dynamical and new 
photometric luminous mass estimates to check the reliability of 
the two methods and to investigate on the IMF of these galaxies. 
Second, we propose to relate tw o samples of non-lens (Coma) and 
lens (SLAC^^ see .Bolton et al.i,200 6) early-type galaxies with rel- 
atively similar physical characteristics. 



3 THE INITIAL MASS FUNCTION 

For the early-type galaxies in the sample, axisymmetric two- 
component dynamical models computed by superposition of or- 
bits ( Schwarzschild 1979) were developed in T07 and T09 to 
reconstruct the galaxy luminous and dark matter composition. 
First, assuming three different inclination angles, each observed 



2D luminosity profile was deprojected to derive a 3D luminos- 
ity profile and convert it into a luminous mass density distribu- 
tion through a constant stellar mass-to-light ratio (T,). Then, a to- 
tal mass density distribution was calculated by adding to the lu- 
minous component a dark matter one, w hich was paramet rised as 
a NF W jNavarro, Frenk & White|[l996l) or a logarithmic l lBinnevI 
Il98lh profile. Finally, the best-fitting model was determined 
following the prescriptions of the maximum entropy technique 
jRichstone & Tremaine 1988.) in superposing orbits in the gravi- 
tational potential associated to the total mass density distribution 
(for further details, see T07 and references therein). The galaxy lu- 
minous mass values M. ^y,,, that are shown here in Table [T] were 
measured from the multiplication of the i?-band luminosities and 
the best-fitting T,, where the minimum-^- value of the stellar mass- 
to-light ratio was chosen between the two different families of dark 
matter haloes. 

We measure the same physical quantity (luminous mass) by 
fitting the galaxy SEDs with composite stellar population mod- 
els. The SEDs are retrieved from the SDSS archive and are com- 
posed of 5 bands extending from 354 to 913 nm. The extinction- 
corrected total magnitudes (modelMag) correspond to magnitudes 
measured through equivalent apertures in all the bands and thus 
provide unbiased galaxy colours, in the absence of colour gradi- 
ents. We u se Maraston's COOSi : indexed M) and Bruzual & Char- 
lot's ( I2OO3I : ind exed B C) dust-free templates a t solar metallicity and 
adopt Salpeter ('1955'; indexed Sal), Kroupa l l200ll: indexed Kro), 
and Chabrier (2003, indexed Cha) IMFs. We assume a delayed ex- 
ponential star formation history (SFH) and perform the SED fits 
on a three-parameter (age, star-formation time-scale, and luminous 
mass) grid of models (for further details, see G09a). The best- 
fitting luminous mass values for Maraston's templates and Salpeter 
(Mf '•") and Kroupa (Mp °'^) IMFs are shown in Tabled 

In Fig.[T] we plot for the 13 galaxies of our sample the dynam- 
ical and photometric luminous mass estimates reported in Table[T] 
In Table |2] we show the average value of the dynamical over pho- 
tometric luminous mass ratio (M,,dyn/W.,phot) obtained with all the 
different templates and IMFs. We observe from Fig. [T] and Tables 
[T] and [2] that the models with a Kroupa and Chabrier IMF give lu- 
minous mass measurements that are significantly smaller than the 
dynamical ones. In particular, the average ratio of the latter over 
the former deviates from one by more than 5cr. On the contrary, 
the models with a Salpeter IMF provide luminous mass estimates 
that are in general consistent with those determined from dynam- 
ics. Despite the overall good agreement of these two luminous mass 
diagnostics when a Salpeter IMF is adopted, we note that for some 
galaxies in the sample the dynamical mass values are moderately 
larger than the photometric ones. 

Although the SDSS photometry is very accurate, the large spa- 
tial extension of the light distribution of the galaxies in the sample 
and the assumption of a de Vaucouleurs profile to estimate the total 
magnitudes may r equire significant corr ections in the photomet- 
ric quantities (e.g., lKormendv et al.ll2009l) . By comparing our best- 
fitting values of the S-band l uminosity with tho se taken from the 
public database HyperLed^ jPaturel et al.l l2003). we observe that 
not all the light belonging to the galaxies may have been included 
in the SDSS photometry. In the simplified hypothesis that the frac- 
tion of light missed in the five SDSS bands is constant and equal to 
that estimated from this comparison in the i?-band, we can correct 
our photometric luminous mass estimates. This results in the empty 
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Figure 1. Comparison of the luminous mass estimates detemiined from 
two-component dynamical modelling (M, jy„) and SED fitting (M» phot) for 
a sample of 13 early-type galaxies of the Coma cluster The values mea- 
sured using Maraston's dust-free templates at solar metallicity and assum- 
ing a Salpeter and a Kroupa IMF are represented, respectively, with squares 
and diamonds. The empty squares represent the Salpeter photometric lumi- 
nous mass values after a simple luminosity correction performed with the 
HyperLeda database. The lo" error bars and the one-to-one congelation line 
are shown. 

squares plotted in Fig.[T]and in average values of M,,dyn/A^*.phot of 
1.14 ± 0.10 and 1.67 ± 0.15 for Salpeter and Ki'oupa IMFs, respec- 
tively, adopting Maraston's templates. 

A reason why the luminous mass measurements obtained 
through dynamical modelling can be slightly overestimated is the 
fact that a small fraction of mass of the dark matter component may 
have been associated to the luminous one. In fact, for a given set of 
observations of a galaxy, if the density distribution of its dark mat- 
ter component resembles closely that of its luminous component, 
then the dynamical models may fit equally well by overestimating 
the luminous mass term at the expense of the dark matter one. Re- 
lated to the previous point, the assumption on the parametrisation 
of the dark matter density distribution have a small, but not neg- 
ligible, influence on the luminous mass estimates. For instance, in 
T07 and T09 the dynamical models with a logarithmic dark mat- 
ter halo provide best-fitting values of the stellar mass-to-light ratio 
that are always equal or larger than those with a NFW halo and the 
majority of the dynamical luminous mass of the sample galaxies 
were derived from the logarithmic best-fitting Y,. On the opposite 
side, possible underestimates of the photometric luminous mass de- 
terminations may be related to dust extinction or metallicity values 
lower than the solar one in the galaxies of the sample. These two 
effects would result in lower IR fluxes tha t in turn would give lower 
mass measurements. Several studies (e.g. jRettura et alj2006h have 
however tested the validity of the dust-free and solar-metallicity ap- 
proximation. We remark that the discussed limitations in dynamical 
and SED modelling may justify the differences observed by com- 
paring dynamical and photometric mass estimates with a Salpeter 
IMF, but not the much larger discrepancies obtained with a Kroupa 
or Chabrier IMF. 



1.4 



1.2 ^ 




0.4 - 



0.2 L , ^ ^ , ^ , ^ \ ^ , ^ J 

10.5 11.0 11.5 12.0 

log(M.,„) 

Figure 2. The B-band stellar mass-to-light ratio versus the luminous mass 
estimates at redshift z = for the early-type galaxies of the Coma (squares) 
and SLACS (diamonds) samples. The values and their lo" error bars are 
measured by using Maraston's dust-free templates at solar metalhcity and 
assuming a Salpeter IMF. The empty squares represent the photometric lu- 
minous mass values of the Coma galaxies after a simple luminosity correc- 
tion performed with the HyperLeda database. The best-fitting correlation 
lines for the two samples are shown. 

4 THE TILT OF THE FUNDAMENTAL PLANE 

The time evolution of the best-fitting SFHs (given by the param- 
eters of age and star-formation time-scale of the best-fitting mod- 
els for the SEDs) allows us to measure for all the sample galaxies 
their luminous mass (M, o) and stellar mass-to-light ratio in the B- 
band (M».oig q) at redshift z = 0. Since the Coma galaxies are local 
galaxies, the values of these two quantities do not differ signifi- 
cantly from those measured at the actual galaxy redshifts. In Table 
[T] and in Fig. |2l we show the luminous mass and stellar mass-to- 
light ratio values that we obtain with Maraston's templates and a 
Salpeter IMF. 

Although the points are relatively scattered, we observe in Fig. 
[2]that the two quantities are correlated: very massive galaxies dis- 
play larger stellar mass-to-light ratios than less massive galaxies. 
The values reported in Table|2]of the Pearson linear ( p) and Kendal l 
rank (f) correlation coefficients (for definitions, see lSalkin3l200d) 
prove that the luminous mass and the stellar mass-to-light ratio are 
statistically correlated at more than a 99% confidence level (linearly 
at a 95% confidence level). We note that the presence of a larger 
scatter at the low mass end is connected to differences in the galaxy 
stellar ages, as highlighted in T09. In passing, we mention that even 
if less evident the same correlation is found between the luminous 
mass and stellar mass-to-light ratio values determined in T09 from 
dynamical modelling. If we fit the scaling relation M.^oL^g ^tc 
as done in Fig. |2l we obtain that our photometric measurements 
give a value of the exponent 6 of 0.22 ± 0.10 (see Table O. As 
reported in Table |2] and shown in Fig. (2] we remark that this last 
result is not sensitive to the adopted IMF and only minimally in- 
fluenced by the luminosity correction that we have performed with 
the HyperLeda data and discussed in the previous section. 



4 C. Grillo and R. Gobat 

Table 2. Important quantities derived from different composite stellar population models and IMFs. 





Sal,M 


Sal,BC 


Kro,M 


Cha,BC 


<M.,dy„/M..pho,> 


1.34 ±0.13 


1.43 ±0.16 


1.96 ±0.19 


2.33 ± 0.24 


e[log(M..o),log(M..oLg'o)] 


0.53 (0.06) 


0.54 (0.05) 


0.53 (0.06) 


0.57 (0.04) 


f[M.,o,M..oLg'(,] ■ 


0.72(<0.01) 


0.66 (<0.01) 


0.69 (<0.01) 


0.55 (<0.01) 


6 


0.22 ±0.10 


0.16 ±0.07 


0.21 ±0.10 


0.18 ±0.08 



Notes - For the the Pearson linear q and Kendall rank f coirelation coefficients, we show their values and, in parentheses, the probabilities that an equal 
number of measurements of two uncon'elated variables would give values of these coefficients higher than the measured ones. 



Based on 6-band measu rements presented in 
|j0rgensen. Franx & K iaergaardI j 19961) and on the FP analy- 
sis performed by'sender et alj (1 1998 1) in a much larger sample of 
Coma galaxies, Treu et al.l | |2006| ) finds that the effective (dynam- 
ical) mass-to-light ratio increases with the effective (dynamical) 
mass raised to the power of 0.21 ± 0.02. From directly observed 
quantities in a sample of nearly 9000 local early-type galaxies of 
the SDSS. lBernardi et al.l Jiooa) also concludes that the dynamical 
mass-to-light ratio scales as the dynamical mass to the power of 
0.22 + 0.05. This relation is often used to parametrise the tilt of 
the FP. The extremely good agreement on the scaling relations 
between stellar and dynamical mass-to-light ratios and masses 
indicates that stellar population effects are the primary source of 
the tilt of the FP. Thi s conclusion is also supported by the study of 
lAllansonet^ ( l2009h on the FP of Coma galaxies. 

By looking at the values determined by T09 of the fraction of 
dark matter over total mass and total mass included within a sphere 
of radius equal to the effective radius of our galaxies, we do not 
find a statistically significant correlation. Hence, an increment of 
the dark matter fraction going from massive to very massive galax- 
ies does not seem to play an important role in defining the tilt of 
the FP. Just a secondary role may also be associated to a deviation 
from homology in the structural properties of early-type galaxies, 
as recently demonstrated by two s tudies on diff erent samples of the 
SAURON dCappellari et"ai1l2006l) and SLACS jSolton et alj2008l) 
surveys. 



5 COMA AND SLACS GALAXIES 

The galaxies of the Coma sample studied here result to be compa- 
rable to those of the SLACS sample analysed in G09a in terms of 
luminous mass properties and as far as the main origin of the tilt of 

the FP is concerned. 

In agreement with the results o f iGrillo et al.l ( l2008j) and G09a 
obtained from two samples of luminous and red lens early-type 
galaxies with luminous mass values in the same range of our Coma 
galaxies, we find that the choice of Maraston's or Bruzual & Char- 
lot's templates does not change remarkably the photometric mass 
measurements, as long as the same IMF is adopted. This means 
that systematic differences between these two models due to, for 
example, their different treatment of the late stages of stellar evo- 
lution do not play a significant role for these samples. Moreover, 
our photometric luminous mass measurements based on the accu- 
rate SDSS multi-colour photometry are consistent, within the er- 
rors, with those obtained from Schwarzschild and stellar dynamics 
plus lensing modelling only if a Salpeter IMF is chosen. We men- 
tion that for old galaxies (like those in our samples) a Salpeter IMF 



and the "bottom-light" IMF proposed bv lvanPokkun] j2008t) are 
in good accord on the values of the stellar mass-to-light ratio. 

After correcting for passive evolution, 15 lens galaxies of the 
SLACS survey were proved to occupy the s ame FP as non-lens 
galaxies of the Coma cluster l lTreu et al.l200 6). For a bigger galaxy 
sample from the same lensing survey, with luminous mass values 
in the same range of the Coma galaxies studied here (see G09a), we 
show in Fig.|2]the physical quantities already plotted there for the 
Coma galaxies. The two samples reveal almost identical scaling re- 
lations between stellar mass-to-light ratio and luminous mass, sug- 
gesting a common explanation for the tilt of the FP. The interesting 
offset between the Coma and SLACS early-type galaxies may be at 
least in part associated to the different environmental conditions in 
which the galaxies assembled their stars. In fact, the SLACS galax- 
ies are mainly field galaxies as opposed to the cluster galaxies of 
Coma. The impact of an older age or a shorter burst of star forma- 
tion in giving higher values of the mass-to-light ratios for the Coma 
with respect to the SLACS galaxies is not addressed here and could 
be checked quantitatively by inspecting the galaxy spectra. 



6 COMMENTS AND CONCLUSIONS 

We have selected a sample of 13 bright Coma galaxies, object of a 
previous detailed dynamical study to differentiate between the lu- 
minous and dark matter components. We have modelled the galaxy 
SEDs (from the SDSS) with different dust-free solar-metallicity 
composite stellar population models and derived luminous mass 
estimates as a function of the assumed IMF. By following the time 
evolution of the best-fitting SFHs, we have also measured the val- 
ues of the galaxy luminous mass and fi-band stellar mass-to-light 
ratio at redshift z = 0. We have concluded that: 

- Given the errors, a Salpeter IMF provides luminous mass es- 
timates that are consistent with the dynamical ones, as opposed to 
the approximately two times smaller values obtained by choosing a 
"top-heavy" IMF such as Kroupa or Chabrier IMFs. 

- The luminous stellar mass-to-light ratio scales with the lumi- 
nous mass like the corresponding effective quantities are observed 
to do, thus suggesting that the tilt of the FP can be mainly ascribed 
to stellar population effects. 

A comparison of the photometric and FP analyses per- 
formed on this sample of Coma galaxies and different samples 
of lens galaxies of the SLACS survey (G09a; Grillo et al. 2008^; 
iTreu et '^ l2006h seems to confirm that old, massive early-type 
galaxies prefer a Salpeter IMF to either a Kroupa or Chabrier IMF 
and exhibit stellar population properties that are responsible for the 
the most part for the tilt of the FP. We remark that early-type galax- 
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ies that are different in age and mass from the galaxies of our sam- 
ple may deviate from the picture presented in this analysis. 
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